ODD PERFECT NUMBERS HAVE AT LEAST NINE DISTINCT 

PRIME FACTORS 



PACE P. NIELSEN 

Abstract. An odd perfect number, N, is shown to have at least nine distinct 
prime factors. If 3 | Af then N must have at least twelve distinct prime 
divisors. The proof ultimately avoids previous computational results for odd 
perfect numbers. 



1. Introduction 

A perfect number is one where (j{N) = 2N. In other words, the sum of the divi- 
sors of N is twice N. These numbers have been studied since antiquity. It is known 
that N is an even perfect number if and only if = (2^ — 1)2^"^ with 2^ — 1 prime. 
Sufficiency was proven by Euchd and necessity by Euler. A prime number of the 
form 2^ — 1 is called a Mersenne prime, and there are currently 43 known. There is an 
ongoing, online, distributed search for such primes at http://www.mersenne.org. 

The search for odd perfect numbers has not been as successful. Currently there 
are none known, making their existence the oldest unanswered question in mathe- 
matics. However there are a great number of necessary conditions, which go through 
periodic improvements. The list of conditions given here is the same list as in |25| . 
but with recent improvements included. 

Let N be an odd perfect number (if such exists). Write N — Y[i=iPT where 
each Pi is prime, pi < p2 < ■ ■ ■ < Pk, and k — lu{N) is the number of distinct prime 
factors. The factors are called the prime components of N. Then: 

• Eulerian Form: We have N = 7r"m^ for some integers tt, Q!,to G Z_|_, with 
TT = a = 1 (mod 4) and tt prime. The prime tt is called the special prime 
of A^. 

• Lower Bound: Brent, Cohen, and te Riele |2| using a computer search 
found that N > 10^°". William Lipp, using the same techniques, is close 
to pushing the bound to > 10^*^°, and plans to start a distributed search 
at the website http : //www . oddperf ect . org. 

• Upper Bound: Dickson [7] proved that there are finitely many odd perfect 
numbers with a fixed number of distinct prime factors. Pomerance |24) 
gave an effective bound in terms of k. This was improved in succession by 
Heath-Brown [El, Cook 0, and finally Nielsen P| to 2^^" . 

• Large Factors: Jenkins ^Hl proved that pk > 10^, and lannucci jl3 | .jl4 | 
proved > 10^ and > 10^. 
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• Small Factors: The smallest prime factor satisfies pi < |A: + 2 as proved 

by Griin g^. For 2 ^ i ^ 6, Kishore [HI showed that pi < 2'^'''\k + 
and this has been slightly improved by Cohen and Sorli 5 . 

• Number of Total Prime Factors: Hare proved that the total number 
of (not necessarily distinct) prime factors of N must be at least 47. In 
unpublished work he has improved this to 75. 

• Number of Distinct Prime Factors: Chein ^ and Hagis (21 independently 
proved that uj{N) ^ 8. Hagis [HI and Kishore [IHI showed that if 3 >; N 
then Lu{N) ^11. This paper improves both of these bounds by 1. 

• The Exponents: For the non-special primes, pi, write Oi = 2bi. If d = 
gcd(26,; + 1) then d^O (mod 3) by a result of McDaniel [T^ . 

With such a number of conditions, it might seem that an odd perfect number 
could not exist. Pomerance has given an interesting heuristic, available at Lipp's 
website, suggesting that odd perfect numbers are very unlikely. 

I would like to thank William Lipp for providing some of the factorizations given 
in the lemmas. Also, much of the terminology, notation, and lemmas of the early 
sections of this paper match those found in |25j (which in turn match those in , 
and earlier work like in an effort to establish a sense of both continuity and 
improvement, and Fd like to thank John Voight for some interesting conversations 
on topics related to his paper. 

2. Fixed Notations and Conventions 

Let N be an odd perfect number and k — lo^N). We will write n for the special 
prime. By N we mean the non-negative integers and by Z+ the positive integers. We 
will use an algorithm to prove our results. At each stage in the algorithm there will 
be prime divisors of N that are known and some that are unknown, meaning that the 
prime divisors are either specified by the algorithm or they are not, respectively. 
This set of known primes will change at every stage of the algorithm as it runs 
through different cases, and so the known and unknown primes are constantly 
changing. A more formal definition will be given in a later section. We let fci be 
the number of known, distinct prime divisors of N (at any given stage), and let 
^2 = fc — fci be the number of unknown, distinct prime divisors. Among the known 
prime divisors of N, some of the prime components are also known (again, known 
being a technical term meaning specified by the algorithm). In other words, if p is 
a known prime divisor of N and if our algorithm yields some a G Z+ so that p"! |A^, 
we say p° is a known prime component. We let £i be the number of known prime 
components of N, and let £2 = k — iihe the number of unknown prime components. 

A word of warning: In some theorems we will assume p is a prime with |A^, but 
this doesn't even mean that p is a known prime, let alone that the prime component 
is known. Further, even if p is a known prime, the component may still be unknown 
(in this technical sense) even though we (from our position) are given that p'^jjA^. 
In other words, we are given the hypothesis p'^WN but the algorithm may not have 
specified p or a. We will, throughout, only use the words "known" and "unknown" 
in the sense of known to us through our algorithm, rather than by hypothesis. 

3. Cyclotomic Integers 

The equation a{N) = 2N can be (trivially) rewritten as a{N)/N = 2 which tells 
us that each odd prime divisor of a{N) must somehow divide N, and vice versa. 
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Thus, we want to study the prime factorization of 

i=l i=l 

Letting $„(x) be the nth cyclotomic polynomial (i.e. the minimal polynomial over 
Q for a primitive nth root of unity), we have the partial factorization 

p"-i = []a>rf(p) 

d\n 

and so 

(1) a{p-')=p—^= n Mp)- 

d\n,d>l 

We are further interested in the factorization of If c and d are integers with 

d > I and gcd(c, c?) = 1 we write Od{c) for the multiplicative order of c modulo 
d. If p is prime, we write Vp for the valuation associated to p. In other words, for 
n e Z_(- we have The following results of Nagell are fundamental and 

are often left as exercises in modern abstract algebra books. 

Lemma 1 f |25[ Lemma 3]). Let m > 1 be an integer, and let q be prime. Write 
m = q^n with gcd{q,n) = 1. 
7/6 = then 

^rn{x) = (mod q) 

is solvable if and only if q = 1 (mod m). The solutions are those x with Oq(x) = m. 
Further, Vq{^m{x)) = Vq{x"^ — 1) for such solutions. 
Ifb^O then 

^,n{x) = (mod q) 
is solvable if and only if q= \ (mod n). The solutions are those x with Oq{x) — n. 
Further, if m > 2 then Vq{^rn{x)) — 1 for such solutions. 

We have immediately from Lemma ^a-nd Equation ^ 

Lemma 2 ( \1'6[ Equation 4].|25l Lemma 4]). Let p and q be primes, q ^ 3, and 
a G Z+. Then 

(vq{p°^^P'> -1) + Vq{a + 1) if Oq{p)\{a + 1) andoq{p)^\, 

Vq{a{p''))^\vq{a+l) lfOq{p)^l, 

I otherwise. 

It turns out that the first case of Lemma (when 6 = 0) is the more common 
means of obtaining factors for odd perfect numbers. However, to make sure we can 
always reduce to that case we need the following result, usually attributed to Bang 
but given other proofs such as in 

Lemma 3. Let m,x E TLj^ with x ^ 2. Then $m(x) is divisible by a prime q with 
0q{x) = m, except if x — 2 and m = 1 or 6, or if x — 2^ — 1 (for some i G Z_|_ ) and 
m = 2. 

Note that we are only interested in this lemma when x = p is a prime dividing 
an odd perfect number N. Thus, the case x = 2 never happens. Also, if to = 2 
this corresponds to the special prime (in the Eulerian form) so a; = tt = 1 (mod 4) 
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and hence cannot be of the form 2^ — 1. Thus, both exceptions in the lemma do 
not affect our work. 

Define (Ti{n) = J2d\n'^^' ^'^^ i 6 Z and n e Z+. It is clear that each of these 
functions is multiplicative, cti = cr is the usual sum of divisors function, and ctq is 
the number of divisors function. The following is immediate: 

Lemma 4. Let N be an odd perfect number. //p"||iV, where p is prime, then for 
each d\{a + 1) the number is divisible by a prime q with Oq{p) = d and q = I 

(mod d). In particular, <j{p'^) has at least ao{a + 1) — 1 distinct prime factors in 
common with N . 



A prime q is called a Fermat prime if it is of the form q — 2^ + 1 for some 
j G Z+. One can show it is necessary that j — 2* for some i G N. It is easily seen 
that if i = 0,1,2,3,4 then 2^' + 1 is prime (i.e. q = 3,5,17,257,65537), but no 
other Fermat primes are known. These primes play a special role in the study of 
odd perfect numbers. This is because the prime factorization of q — 1 is exactly a 
power of 2. 

The first thing we can do is restate Lemma in terms of divisors of iV, and 
Fermat primes. 

Lemma 5 f |23|.|25l Lemma 5]). Let N be an odd perfect number, p°'\\N with p 
prime, and let q be a Fermat prime. Then: 

( Vg{p + 1) + Vg{a + 1) ifTT=p=-l (mod g). 



Proof. First suppose that p = 1 (mod q). In this case Oq{p) = 1 and so we just use 
Lemma |21 

Next suppose p ^ 1 (mod q). Note Oq{p) must be a divisor of g — 1 = 2^ . 
However, from the Eulerian form, we see that a + 1 is divisible by 2 if and only if 
p = TT. Further, in this case a = 1 (mod 4) and hence 2||(a + l). Thus Vq{a{p'^)) = 
unless p — t: and Oq{p) = 2. But this last equality is equivalent to p = —1 (mod q). 
Lemma El gives us the needed valuation in this case. □ 

Lemma |31 tells us that there are only a few sources in cr{N) for copies of q. In 
particular, if q^\N for some large n, we can force the size of the special prime to 
be large. To prove this we first need another lemma. 

Lemma 6 ( \25\ Proposition 9]). Let N be an odd perfect number, let q be a Fermat 
prime, and suppose p°'\\N with p prime. 

(i) Ifp 7^ vr and q^\a{p"') then a {p"") is divisible by b distinct primes ri,r2, . ■ .,ri, 
with ri = 1 (mod q*). 

(ii) If p = TT, p = —1 (mod q), and q'^Ka + 1), then <j(jp°') is divisible by 2c 
distinct primes ri,r[, . . . ,rc,r'^ with ri = r[ = 1 (mod g*). 

Proof. For part (i), by Lcmma|3wc have q^\{a + l). So we take to be the divisor 
of ^qi{p) specified by Lemma 0] Part (ii) follows from the same lemmas, noticing 
that since p = tt is special we have 2|(a + 1), and hence we can take and r'^ to be 
the factors specified by Lemma 0] of $gi(p) and ^2q'{p) respectively. □ 



4. Fermat Primes 




ifP=^ (mod q), 
otherwise. 
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Proposition 7 (c.f. |25l Proposition 10]). Let N be an odd perfect number, and 
let q be a Fermat prime with q"'\N. Suppose k, ki, k2, £i, and £2 have their usual 
meanings. Further suppose q^\\(j [known prime components of N). Finally, let k[ 
(respectively £[) be the number of distinct prime factors among the fci known prime 
divisors (respectively, among the £1 knoivn prime components) which are congruent 
to 1 (mod q). 
If 

T = n-b- {k[-£[ + k2)(k[ + /C2 - 1) > 

then TT = — 1 (mod q) and tt is among the unknown prime components. If, further, 
each known prime, p, with unknown component and with p = I (mod 4) satisfies 
p ^ —1 (mod q'^), then tt is among the unknown primes, and 

VqiiT + l)^n-b-{k[~ e[ + k2~ l){k[ +k2-2)- l{k[ + k2- 1)/2J = t'. 

Proof. First, we will show the contrapositive of the initial statement. Suppose tt is 
among the known components, or if it isn't then tt ^ — 1 (mod q). Let p"" be an 
unknown component of N. Lemma |S1 combined with what we've just said, implies 
Vq{a{p"')) = Vq{a + 1) ifp = 1 (mod q), and equals otherwise. Using the equation 
a{N)/N = 2, there are n — b copies of q that must be accounted for by a of the 
unknown components. Note that there are at most k[ + ^2 distinct prime factors of 
N that are congruent to 1 (mod q). Thus, by LemmaEl at most k[+k2 — l copies of 
q divide a+ 1 if p = 1 (mod q) (since p itself cannot divide cr{p°-)), otherwise we end 
up with too many distinct prime divisors of N which are congruent to 1 (mod q). 
There are at most k[ — £{ + k2 primes p with unknown component satisfying p = 1 
(mod q). Hence at most {k[ — £\ + k2){k!^ + ^2 — 1) copies of q can be accounted 
for by a of the unknown components. Therefore 

T^n-b-(k{-£!^ + k2){k[ + fc2 - 1) 

or we will have left over copies of q not accounted for in cr[N). 

To get the last statement, now suppose r > 0, tt = — 1 (mod q), and tt is among 
the unknown components. Lemma |S1 tells us that the only place the extra r copies 
of q can be accounted for is tig(7r + 1), and hence q^Kvr + 1). Therefore tt = — 1 
(mod q'^) and so tt is an unknown prime. 

This means that when we counted the maximum number of possible (unknown) 
prime divisors = 1 (mod q) we included one too many. Thus there are at most 
k[ + fc2 — 1 primes dividing N that are congruent to 1 (mod q), and hence at most 
{k'l — £[ + fc2 — l){k'i + — 2) copies of q can be accounted for by a of the unknown, 
non-special components. For the special component tt", by Lemma |H1 part (ii) we 
see that at most {k[ + k2 — l)/2 copies of q can divide a + 1, else we obtain too 
many factors of TV congruent to 1 (mod q). Since, in fact, an integer number of 
copies of q divides a + 1, we can take the floor. Putting this all together, r' copies 
of q must divide tt + 1 . □ 

Thus, if a large power of a Fermat prime divides N, we see that tt + 1, and 
hence tt, must be large. This isn't quite good enough to simplify our search to 
manageable cases. We need a way to find another large prime divisor of N. The 
trick is to consider divisors of a{q"). Suppose g is a Fermat prime, with g"||7V, and 
n very large. Then using the previous lemma, we can force tt to be very large. One 
might wonder if 7r|cr(g"). The following result speaks to this issue. 
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Lemma 8 ((2l Lemma 1]). If p and q are odd primes with p\(T{q'') and q"^\{p + 1) 
then k ^ 3m. 

So, using the terminology of Proposition [3 if r > 0, g"^ | (p + 1) for the known 
primes, and also 3t' > n, then tt \ cr[q") by Lemma |H1 But since n is big we would 
expect a large prime divisor of N which divides a{q^). The following work clarifies 
how large a divisor we can find for (j{q^). 

Lemma 9. Let p he an odd prime and let q — 3 or 5. If q^^^ = 1 (mod p^) then 
either {q,p) = (3, 11) or q°p^'i^ — 1 has a prime divisor greater than 10^^. If q = 17 
then either {q,p) = (17,3) or q°p^i^ — 1 has a prime divisor greater than 10^^. 

Proof. The papers "^tF and TSJ give a list of {q,p) for which qP^^ = 1 (mod p'^) 
and p < 10^^ with g = 3 or 5 (or p < 10^^ with q = 17). In the cases {q,p) ^ 
(3, 11), (17, 3) the following table gives the requisite factor of g°p(9) — 1. 



9 


V 


Large factor of g°p('?) — 1 


3 


1006003 


154680726732318637 


5 


20771 


625552508473588471 


5 


40487 


625552508473588471 


5 


53471161 


50493456782731 


5 


164533507 


52082118058261 


5 


6692367337 


8930008316757509 


5 


188748146801 


40093613041379 


17 


46021 


1365581260423071390161 


17 


48947 


63895279579889 



□ 

For use shortly, we make the following definition. Letting p, q be odd primes, 
p ^ q,^e set 

if2to,(p), 

or 4 1 Oq {p) and either 
Oq{p) (i)p = 7ror 

(ii) p = 1 (mod 4) and tt is not 
among the known components 
otherwise. 



In other words, o'q{p) is the usual order function, unless it is impossible for both 
p°-\\N and Oq{p) \ (a + 1) to hold, due to consideration of the Eulerian form. 

Proposition 10. Let N be an odd perfect number with k, ki, and k2 having their 
usual meanings. Suppose q — 3 or 5 is a known prime divisor of N , q^^\\N , q ^ ■n 
and TT I (t(<z"). Suppose pi, . . ■ ,Pki-i are the other known factors of N, besides q. 
For each i e [1, fci — 1] define 

'O lfo'p^{q) = 

= { max(,s + < - 1, 1) if o'p^{q) ^ 0, s ^ Vp,{a{q°p^^^^~^)), 

and t e Z+ minimal so that p\ > 100. 

Set V = YiiLi ^ Vi ■ Suppose tt is among the fc2 unknown prime factors, and ^2 > 1. 
Finally, assume that all unknown prime factors are greater than 100. 
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II 



V \ V 



then (T{q"') has a prime divisor among the unknown primes at least as big as the 
above minimum. If q ~ 17 one can replace 10^^ with 10^^, and then the result still 
holds. 

Proof. Wc do the case when q ~ 3, since the other cases are similar. First suppose 
(T(g") = ((7"+^ — l)/{q — 1) is at most divisible by p,^' for the known primes, and 
square-free for the unknown primes. Then since tt, g | cr(g"), the largest unknown 
divisor of (T(g") is at least 

'a{q'')\T^ 



V 

unless this quantity is ^ 1 (in which case there might be no unknown factors). 

So we may assume there is some prime p\N, o'p{q) ^ 0, so that (T{q^) is divisible 
by lip is unknown, or p*^"*"^ if p is known (and e is the corresponding ti), with p 
maximal among such primes. By Lemma^lwe may also assume that if \ (g^^^ — 1) 
and p is unknown then p|(n + 1). (This is where 10^^ comes into the minimum.) 

Thus in either case, p*|(n + 1), where > 100 (taking t = 1 if p is an unknown 
prime). Then we have 

a(/-i)|a(g"). 

Thus it suffices to find a large divisor of {(f^ — l)/(g — 1). By Lemma ^ (jf^ ~ 
l)/{q— 1) is only divisible by primes larger than p, or p itself to the first power. 
(In this case, q being Fermat means the quantity isn't divisible by p, and we could 
replace max(s + t — 1, 1) by s + t — 1 in the definition of e^. But to keep similar 
notations later when we take q to be an arbitrary prime, we don't use this fact.) 
But then, by the maximality condition on p, [qP — l)/(g — 1) is not divisible by 
more than pj' for known primes and the first power for all the unknown primes. So 
the analysis we used in the first paragraph goes through by only changing n + 1 to 
p*. Finally, note that p* > 100, so we have the appropriate bound. □ 

The most useful case when we will use Proposition llOl is when n is very large and 
ki is close to k. So, in practice, we will usually end up with 10^^ as the lower bound 
on a divisor of ct((7"). In fact, a lot of work could be saved if the bounds given in 
[TB| were improved, or there was some means to work around square divisors. 



5. Non-Fermat Primes 

Sometimes our odd perfect number will not be divisible by a large power of a 
Fermat prime, but rather by a large power of some arbitrary prime, q. Unfortu- 
nately, we don't have Lemma and so we can't put all of the "extra" factors of q 
into 77 -I- 1. This causes two problems. First, we have to spread the extra factors of 
q among the unknown primes, thus reducing the number of extra factors we have 
at hand, exponentially. Second, if p ^ 1 (mod q) is one of our unknown primes, we 
cannot reduce to the case q\'^2{p) but rather q\'^d{p) for some (arbitrary) d > 1, 
d\{q — 1). Thus, we need a way of bounding the size of p for which (7" |$d(p)- This 
second problem is easily dealt with. 
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Lemma 11 Theorem 2]). Let q be an odd prime. Let ai be a primitive root 
modulo q. For r ^ 2, define — a\ (mod q'). Then {a™ (mod q^) \ m — 
0, . . . , g — 2} gives a complete set of incongruent solutions to a'^^^ = 1 (mod q^). 

Lemma 12. Let q < 1000 be an odd prime. Ifp'^~^ = 1 (mod q") for some n G Z+ 
and some odd prime p, then p ^ min(g"~^, 10^") except when 

{p,q) = (40663372766570611389846294355914421,7). 

Proof. Let q < 1000 be an odd prime. Fix m G N so that ^ j^qSO ^^^j 

is minimal. A computer search using Lemma |l II demonstrates that every positive 
solution to x'^~^ = 1 (mod g™) with x 1 satisfies x > 10^°. Hence the same 
is true if we replace m with a larger integer. Thus, it suffices to search for prime 
solutions to x'^^^ = 1 (mod g") for n ^ m. Again, a computer search yields the 
stated result. □ 

The choice of 1000 in Lemma [TB is easily improved, but is large enough for our 
needs. Also, the exceptional {p, q) in the lemma is irrelevant to our work, since 
in this case we find Oq(ji) — 6 and a(j)^) gives rise to at least 12 additional prime 
factors of N besides p and q. 

Lemma 13. Let N be an odd perfect number, let q be an odd prime, and suppose 
p°'\\N with p prime. 

(i) If p = \ (mod q) and q^\a[p°') then (j{p°') is divisible by b distinct primes 
ri, r2, . . . , with ri = 1 (mod q''). 

(ii) If p ^ 1 (mod q), Oq{p)\{a + 1), and q'^\(a + 1), then a{p°') is divisible by 
c ■ cfid{oq{p)) distinct primes divisors rij, i e [l,c], j G [1, ao{oq{p))], with 



rij — 1 (mod q*). 
Proof. Analogous to Lemma 



□ 



Proposition 14. Let N be an odd perfect number, and let q < 1000 be a prime divi- 
sor of N with q"\N. Suppose b, k, ki, k2, ii, £2, k[, andi'^, have the same meanings 
as in Proposition^ Suppose further that the exceptional case of Lemma \lS\ doesn't 
hold. Let T be the set of known primes with unknown component, different from q, 
and ^ 1 (mod q). Let 

k[ + k2 



per, o'(p)^o 



(^q{p°-^P^ - 1) + 



cro(og(p)) 



-(fci-€'l+fc2)(fci+fc2-l). 

If T > Q then one of the unknown primes is not congruent to 1 (mod q). Further, 



in this case, one of the unknown primes is at least as large as TD.\n[q'^ ^, 10^'^) where 



T — mm 



n — b — 

^ peT,o'(p)#o 



Vqip' 



Oqip) 



-{k[ ~ e[ + k2- m){k[ +k2-m~l) 



1) 



k[ 



k[+k2 



MOq{p)) 



Proof. The proof is similar to Proposition From the equation a{N)/N = 2 we 
know that q"'\a{N), and so we try to account for as many copies of q in a{N) as we 
can. The quantity r is exactly how many copies of q are unaccounted for, if all of 
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the unknown primes arc = 1 (mod q), and we try to account for as many copies of 
q as possible from the known primes using Lemmas El and 1131 If t > this means 
we actually have left over copies of q, which yields a contradiction, and hence not 
all of the unknown primes are = 1 (mod q). 

In this case, let m be the number of unknown primes ^ 1 (mod q). Lemma |5] 
tells us that there are two sources of copies of q in <j{N); the exponents of the 
primes, and p°i'^P'> — 1. The quantity 

~ {k\ - + fc2 - m)(k'^ + fc2 - m - 1) - m + ^2 " m 

is the number of copies of q in (j(iV) not yet accounted for, after we account for 
as many copies of q as we can (again using Lemmas 1^ and [T^ except those copies 
of q which come from cr(p°«^P-'~i) for the unknown primes p with p ^ 1 (mod q). 
Thus, if we divide Equation El by m, and take the ceiling, we have a number of 
copies of q that must be accounted for by p°9(p) — 1 for an unknown prime p with 
Og(p) ^ 1. Since m is unspecified by the hypotheses, we take the minimum over all 
possibilities. Finally, we apply Lemma E| □ 

Next we want to prove a result analogous to Proposition llUI except for arbitrary 
primes. However, there have to be a few differences. First, we need something to 
play the role of Lemma |S| It turns out that it doesn't hurt much to just assume 
that the large prime, call it p, coming from ProDOsition ll4l mav in fact divide cr(g"). 
So we can write (7"+^ — 1 = (q — X)]fm for some c G N and m G Z+. Powering this 
equation to the (g — l)st power, since p'^^^ = 1 (mod (f ) and r' ^ n, we have 

{{q - l)m)«-i = ((g - l)mpy-^ = (9"+^ - l)""! ee 1 (mod 9^') 

and hence we can bound (g — l)m using an analogue of Lemma [T^ (where we look 
for solutions to the equation a;'"^ = 1 (mod 5") which are divisible hy q—1, rather 
than prime solutions). The following lemma provides this result: 

Lemma 15. Let q < 1000 be an odd prime. Suppose a''~^ = 1 (mod g") for some 
n G Z+ and some integer a with (q — l)|a. If q ^ 11 then a ^ min(g"^^, 10^°). // 
q^7, then a ^ iiim{q''-^, 10^°). 

Proof. An easy computer search as in Lemma [T^ □ 

Next we need to work with possible square factors of cr(g"), similar to LemmaEl 

Lemma 16. Let p and q be primes with p G (lO'^, 10^^) and q — 7, 11 or 13. // 
qP^^ = 1 (mod p^) then a{q°''^'^'^^^) is divisible by two primes greater than 10^^. 

Proof. By 16 , there are only 3 pairs (p, q) satisfying the conditions, namely 

(p, q) = (491531, 7), (863, 13), (1747591, 13). 

In the first case, since 65|op(g), just factor cr(7^4) to find two distinct primes greater 
than 10^1. In the second case Op{q) — 862. One factor, 16002623839393, is easily 
found and another is 

14-8633 -68099" 
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In the last case, since 195\op{q), one factors o'(13'^^) to find one prime greater than 
10^^ and factors (t(13^^) to find the other one. (Or, to save time, use the onhne 
factorizations on Richard Brent's extensions of the Cunningham tables for these 
last two.) □ 

We are now ready to prove: 

Proposition 17. Let N be an odd perfect number and let q — 11 or 13 be a known 
prime divisor of N , with q^^\\N. Let t,t' be as in Proposition^^ suppose all the 
hypotheses of that proposition are met, and let p be the guaranteed unknown prime. 
Let pi, ■ ■ ■ ,pki-i be the known primes different from q. Let be defined as before, 
and put V — Yli^i^ pI^ ■ Suppose k2 > 1- Finally, assume that all unknown prime 
factors are greater than 100. 
If 

-l50 \ / „r'-2 \ ''-2-1 



then cr(g") has a prime divisor, different from p, among the unknown primes, at 
least as big as the above minimum. Lf q = 7 the same result holds if we replace 
g^'"^ by g^'"^. 

Proof. We do the case g = 11 or 13, the other being similar. So assume the above 
minimum is > 1. First note that if d\{n + 1) then (j(q'^^^)\a{q''^) and so it suffices 
to show that a{q'^^^) has a prime divisor larger than the above minimum, different 
from p, for some d\{n + 1). 

By the method of proof given in Proposition 1101 we may assume that at most 
Ei copies oi Pi divide (T(g'*~^), for some d either greater than 100 or equal to n + 1. 
Further, because 10^^ occurs in the above minimum, we may assume that the only 
unknown prime that may divide a{q"') which may be greater than 10^^ is p. Then 
by Lemma Il6l and the fact that the unknown primes are greater than 100, we may 
assume (T{q"') is square-free for unknown primes, except possibly p. 

From the proof for Proposition ^1 we have = 1 (mod q'^ ). Write q"^ — 1 ^ 
(q — l)mp'^ with c G N, m G and gcd(p, m) = 1. Powering this equation to the 
{q — l)st power, we have 

{{q ~ l)m)«-l = {{q - l)mpy-^ = [q'^ - If-^ = 1 (mod g«""(^Moo)), 

By Lemma [TCI {q - l)m ^ min(g^'-^ q^^ , 10^°) = min{q'''-'^, 10^"). Thus 

— ^ mm 



\{q-l)V' {q-l)V )- 

Since m/V is at least as big as the part of cr{q'^^^) made up from the unknown 
primes, different from p, if we take the k2 — 1 root of the minimum we have the 
appropriate lower bound. □ 

6. Abundance and Deficiency 

Let n G Z+. Recall the multiplicative function cr_i(n) = X^din '^^^ introduced 
earlier. This function can alternatively be written using the formula a_i(n) = 
a{n)/n, and so a{n)/n = 2 if and only if (T_i(n) = 2. A number n is called 
abundant when a-i{n) > 2 and deficient when a-i{n) < 2. We can use abundance 
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and deficiency computations to limit choices on possible prime factors of an odd 
perfect number N. First, we extend the definition of cr-i by setting 

a_i(p°°)= hm a^,{p^)^^. 



Lemma 18 ( \25i Proposition 2]). Let p and q be odd primes. Ifl^a<b^oo 
then 1 < (T_i(p°) < cr_i(p''). // 1 ^ a, 6 ^ oo and p < q then tT_i(q^) < a_i{p°-). 



Lemma 19. Let N be an odd perfect number. Suppose Pi, ■ ■ ■ ,Pki cire the known 
prime factors of N , p^^\N, anc 
the smallest unknown prime is 



prime factors of N , p^'\N, and ki < k = uj{N). // H = UiLi^^-iipT) < 2, then 



n 

Pki+l > 7^ 



n 

Proof. We find 

' ki \ 

Pki+i + 1 



2 = a_i(iv) Yla-iipT) ^-i{pk,+i) = n • 



pfcl+1 



where the inequality in the middle follows from Lemma ^| Noting 11 ^ 1 , we 

obtain ^ > \ -\ — . Therefore ^ — - — and taking reciprocals gives us the 

n ^ Pfci+i n ^ Pki+i t, f t> 

result, since 2 — 11 > 0. □ 



Note that in the lemma if 11 > 2 then Y]!iLi pT abundant, hence N is abundant. 
If n = 2, then Hiili pT already an odd perfect number. 

The following lemma is the true key to our search for odd perfect numbers, as 
simple as the proof is (after wading through the hypotheses). This is because we 
built up machinery in the last few sections to find bounds for large prime divisors 
of iV. 

Lemma 20 (c.f. |5| Lemma 2.2]). Let N be an odd perfect number. Let pi, . . .pk 

be the prime divisors of N , and let Oi be such thatp1^\\N. Fix the numbering on the 
indices so thatpi, . . . ,pi^ are the primes with known prime component, pi+i, . . .pki 
are the other known primes, andpk^+i < . . . < Pk are the unknown primes. Suppose 
among the unknown primes we have bounds pk > Pi, . . . ,Pk-v+i > Pv, with v < k2 
and Pi > 1 for each i £ [1, v]. For each u G [0, v], set 



it^-ibD 




Finally, suppose k2 > 0. 

If Au < 2 then the smallest unknown prime is 

Au{k2 -u) 

Therefore, 

' Au{k2 - u) 



Pki+i ^ min 



e[o,i>], A„<2 \ 2 - A, 
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Proof. We compute 

fc 

2 = a_i(iV) = n'^-ibr) 
1=1 

/ il \ / fe-u \ / ^ 

< H'^-i^^'D n '^-i(pr) n 

fc — A;— ii — A,'i— 1 

i=A;i + l i=0 
_ . Pfci+1 + fc-M-fcl-l _ / k2 -U 



Pfci+1 - 1 V Pfei+1-1, 

Now, recall that u ^ v < k2 which implies k2 — u ^ 1. Also < A„ < 2, so we 
solve the main inequality as we did in the previous lemma, finding 

A„(fc2 - u) 

The last statement follows. □ 

One major diff'erence between Lemma 1^ and Lemma [T^ is that if A„ > 2 then 
that doesn't necessarily imply N is abundant. (It is true that if /ci = fc and Aq < 2 
then N is deficient, however.) This means that we might end up with Aq > 2, and 
hence we have no upper bound on Pki+i- 

7. The Algorithm 

To verify that an odd perfect number has at least 9 distinct prime factors, we use 
a factor chain algorithm to check all possible cases for odd perfect numbers with 
exactly 8 distinct prime factors, and find contradictions in each case, similar to the 
algorithms described in both |S] and This section will describe the algorithm, 
and the next section will explain my implementation. 

We start by knowing 3|A^, so we consider each of the cases 3^||A^, S^'IjA^, and so 
forth. We think of each of these cases as branches on a tree. On the branch 3^||A^ 
we know that a{3'^)\2N and so 13|A^. Hence the case 3^||A^ branches further and 
we have to consider each of the subcases 13||A^, 13^||A^, 13*||A^, and so forth. We 
continue finding new factors, and the tree continues branching out. 

Some cases do not result in new primes. To continue the factor chain in these 
cases we use the bounds of Section |S1 to find an interval for the smallest unknown 
prime. Then we consider all the new primes in the given interval. We only stop 
branching when we reach a contradiction (such as having too many factors) . 

As it stands there are still an infinite number of branches. To get around this 
problem, we can combine all the cases p^\\N for large n together into one case. 
More precisely, let i? be a large integer (around the size 10^°) which we fix at 
the beginning of the algorithm. Then, once we reach the branch p", with n large 
enough so that > B, we no longer consider the case p"||A^ but rather we just 
assume p°|iV for some a ^ n, and thus consider all the remaining cases together. 
(In particular, p is still a known prime divisor of iV, but the component is unknown, 
although we are given a lower bound of B on the size of the component.) We then 
label this composite branch by p°° . 
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So, for example, if B = 10^ then since 3"* < i? < 3^ the first level of our tree 
would consist of the branches 3^, 3"*, and 3°°. When we are on a branch with p°° we 
say p is an infinite prime (not to be confused with the infinite primes of algebraic 
number theory). Notice that infinite primes do not provide more primes for the 
factor chain, and so we have to rely on on the interval bounds of Sectional If we set 
B too low then the primes become infinite too quickly and we may have the case 
the there is no upper bound for an interval. (This corresponds to the case when 
Ao > 2 in Lemma |5ni) If we make B large enough this never happens (for a proof 
see 121), and the algorithm will only have to consider a finite number of cases. 

We are now able to expressly define what we mean by known primes and compo- 
nents. Suppose we are on the branch 3°°5'*. Then the known primes are 3, 5, 11, 71 
(the primes 11 and 71 come from a{5*)\2N), and the only known component is 5^. 
In this case we say that 3 and 5 are on while 11 and 71 are off. In other words, the 
on primes are exactly the known primes for which we have started the branching 
process. Note that fci — £i is exactly the number of (known) primes which are 
infinite or off. 

In our example, since 11 is the smallest off prime we continue the branching 
process on this prime, rather than 71. When there are no off primes we use the 
interval bounds to arrive at another prime, as explained earlier. Whenever we reach 
a contradiction, we go to the next available branch. 

The following is a possible, initial print-out in the case fc = 5, i? = 50. (Note: 
When calculating intervals this output only used Aq in Lemma EH) 

3-2 => 13-1 

13"! => 2-1 7-1 
7-2 => 3-1 19-1 

19-00 : 21 < p_5 < 23 N 
Too : 9 < p_4 < 21 

11-00 : 374 < p_5 < 540 

The letter N means that there are no primes in the interval, which is a contra- 
diction. Adding more contradictions, or using the full power of Lemma 1201 can 
further simplify the output. 

8. An Implementation 

There are three main differences between our implementation of this algorithm, 
and the implementations in and |25| : 

First, the bound B is not allowed to increase within the algorithm. Allowing B 
to vary fully automates the algorithm at the expense of unnecessary complexity. 
The number B is fixed at the outset, and only increased manually if needed. 

Second, the use of Lemma 1201 allows for stronger upper bounds on intervals. 
In the terminology of that lemma, we always have v ^ 3. We take Pi = 10*", 
P2 = 10^, and P3 = 10^ unless we already have known prime divisors larger than 
these bounds, or unless these bounds are superseded by the work in Sections 01 
and|S| (So, for example, if the largest known prime is > 10^, and the next largest 
known prime is < 10^, then we could take Pi = 10"^ and P2 = 10^ If large powers 
of small primes divide TV then we can use the work in previous sections to increase 
Pi and P2.) 
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Third, the contradictions used are different. Here is a complete list of the con- 
tradictions in our implementation: 

MT There are too many total factors. 

MS There are too many copies of a single prime with known component. 
S There is an off prime smaller than an on prime coming from interval com- 
putations. 
A The number is abundant. 

D There are k known primes, and Aq < 2, hence N is deficient. 

F The special prime n belongs to a known component, but the hypotheses of 

Proposition [7| hold showing tt must be in an unknown component due to a 

Fermat prime. 

N There are no primes in the interval given by Lemmas 1 1 91 and HUl or all the 
known primes are on and the only primes in the interval are already known. 
SFl There are fc — 1 known primes, and the interval formula gives an upper 
bound of pk < C, but we know from the fact that a large power of a small 
Fermat prime divides N that some unknown prime is larger than C, by 
Proposition 13 

SF2 Similar to SFl except we have a contradiction between the interval formula 
and PropositionllUI 

SNFl Similar to SFl, except we have a contradiction from a small non-Fermat 

prime, using Proposition [T^ 
SNF2 Similar, using Proposition El 

Plint There are k — 1 known primes all smaller than 10^, and the interval formula 

gives an upper bound pk < 10''. 
P2Int There are k — 2 known primes all smaller than 10'*, or fc — 1 known primes 

with the largest > 10*" and none other > 10^; and the interval formula says 

the next prime is < 10'*. 
P3Int There are fc — 3 known primes all smaller than 10^, or fc — 2 known primes 

only one greater than 10'* and none other > 10^, or fc — 1 known primes 

with one larger than 10^ and another larger that 10^ and none other > 10^; 

and the interval formula says the next prime is < 10^. 

Of course, before the algorithm checks for any contradictions it always checks if 
we have an odd perfect number. There were other contradictions we might have 
included. For example, if there are fc known primes, and none of them can be 
the special prime, this means N cannot be an odd perfect number. However, the 
strength of our method lies in the fact that we rarely have exactly fc known primes, 
and in those cases one of the other contradictions will do. 

9. Points of Improvement and the Results 

The algorithm was implemented in Mathematica on a Pentium 4 personal com- 
puter. The factorizations of cr(p°) were carried out using a only a probable primality 
test, good for integers < 10*^. Thus, if a factor of <j{p'^) was larger than this bound, 
Mathematica's primality proving routine was also run. After running the algorithm 
it became clear that certain modifications would speed up the process. First, the 
bound B was too uniform. So, B was replaced by two bounds Bi and B2, where a 
prime p < 1000 became infinite when p"-\N and p°- > Bi, while a prime p > 1000 
became infinite when p°'\N and p"^ > B2- 
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With Bi = 10^° and B2 = 10^^, the algorithm didn't terminate until we reached 
the case 

3°°5°°17°°257°°65537°° 

because B2 wasn't large enough. After increasing it to B2 — 10'^" the algorithm 
finished the case u^N) = 8 without finding an odd perfect number, running for a 
total of less than 3 days, after numerous stops and starts. 

One of the longest cases (besides the one mentioned above) was 3'*7°°5^. This 
is because the bound given in Proposition 1171 depends on taking (^2 — l)st roots. 
So, when k2 — 3 the bound is really around the square-root of where it should be. 
This led to formulating the improved result: 

Proposition 21. Assume all of the hypotheses of Proposition^^ hold. Further, 
suppose that we can bound pkj^+i < P. Finally assume k2 > 2. // 

(replacing lO^^ hy 10" if q ^ 17) th en a{q") has a prime divisor among the un- 
known primes at least as big as the above minimum. 

Proof. One reduces to the case that criq"") is square-free for unknown primes, just 
as in PropositionllOl Then the result is obvious. □ 

Of course, Proposition El is similarly improved. Once these improvements were 
put into place, the case 3'*7°°5^ took less than an hour. The case uj{N) = 8 took 
less than half a day of non-stop running, with less than 800,000 lines of output. 
We thus have: 

Theorem 22. An odd perfect number has at least 9 distinct prime divisors. 

Running the algorithm with Bi = 10^", B2 = lO^'', A: = 11, and forcing 3 ] N, 
also terminated without finding any odd perfect numbers after about one million 
lines of output. Thus: 

Theorem 23. An odd perfect number N with 3 \ N has at least 12 distinct prime 
divisors. 

At this point, the results seemed to depend on a lot of previous computation. 
To avoid this, the algorithm was rerun with the following changes: 

• All uses of the Plint, P2Int, and PBInt contradictions were removed. 

• The results of jTlj, and were not used in interval calculations. 

• It wasn't assumed that 3\N. In fact, Griin's result wasn't even used, but 
rather Lemma EUl to obtain pi ^ k + 1. 

• It wasn't assumed that uj{N) ^ 8. So each k £ [1, 8] was checked. 

With these changes, the only results (outside this paper) which were used were (i) 
some of the easy to prove non-computational lemmas (such as the Eulerian form) 
and (ii) the main result of |lf)j for primes a ^ 17. The algorithm again terminated 
after about one million lines of output, verifying Theorem l22l A similar run verified 
Theorem!^ 
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10. Future Results 

To do the case w^N) — 9 with these techniques, in a reasonable amount of time, 
it is necessary to find bounds for pk-2, or increase the lower bound on Pk-^i- For 
example, consider the case of 

3°°5°°17°°257°°65537°°. 

We can make the powers on these primes so large that we may, for all intents and 
purposes, assume pk is as big as we like (except not quite big enough to use 22 ). 
However, the best we can do with our tools for pk-i is 10^^, which is much smaller 
than the lower bound already given by the interval formula. The interval for pk-i 
has a length of about 10^^, which is much too large to check one prime at a time. 
Further, it is currently computationally unfeasible to push the bound in 16 up to 
1019. 

Another possible line of attack would be to suppose that two primes qi , (72 divide 
N , each to a high power. Then one might be able to prove that ^(g"^) and criq-^^) 
must each have large divisors, different from one another, once ni and 112 are large 
enough. Another alternative would be to show that the square-free part of cr(g") 
becomes large as n does. 
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